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ABSTRACT 


The young binary system HD 350731 is a noteworthy laboratory for study¬ 
ing early-type binaries with similar components. We present here the analy¬ 
sis of differential multi-color photometric and spectroscopic observations for the 
double-lined detached system. Accurate absolute parameters were determined 
from the simultaneous solution of light and radial velocity curves for the first 
time. HD 350731 consists of two B8V-type components having masses and radii 
respectively of M l = 2.91 ± 0.13 M Q , M 2 = 2.80 ±0.14 M 0 , R\ = 2.11 ± 0.05 R 0 
and R 2 = 2.07 ± 0.05 R 0 . The effective temperatures were determined based on 
analysis of disentangled spectra of the components and derived to be 12000 ± 250 
K and 11830 ± 300 K for the primary and secondary components, respectively. 
The measured projected rotational velocities, 69.2±1.5 km s^ 1 for primary and 
70.1±1.7 km s^ 1 for secondary, were found closer to the pseudo-synchronous ve¬ 
locities of the components. Comparison with evolutionary models suggests an age 
of 120±35 Myr. Kinematic analysis of the unevolved binary system HD 350731 
revealed that it belongs to the young thin-disc population of the Galaxy. 

Subject headings: stars: fundamental parameters — stars: binaries: eclipsing — stars: 
individual: HD 350731 - Techniques: photometric, spectroscopic. 



1. Introduction 


Precise knowledge of the absolute parameters of stars (mass, radii, etc.) is a means 
to better understand the structure and evolution of galaxies. Eclipsing binary stars are 
the most important objects for determining these basic parameters of stars. In order to 
derive the mass, radius, temperature and other physical parameters of eclipsing binary 
stars, spectroscopic, photometric -and interferometric data are required. In particular, the 
quality of the observational data plays crucial role in determining the absolute parameters 
of the component stars precisely. Among the eclipsing binary stars, detached double-lined 
spectroscopic binaries are the best sources for acquiring more accurately the main physical 
properties of stars (e.g. Southworth 2013, Lacy et al. 2015). Well measured detached 
eclipsing binaries have been recently listed and studied based on distributions of their 
absolute parameters by Torres et al. (2010) and Eker et al. (2014). 

HD 350731 (BD+20 4323, GSC 01624-00493, V = 9 m .60) was identified as an eclipsing 
binary system with eccentric orbit by Otero et al. (2004). In the AGK3 catalogue, the 
spectral type of the system is given as A0 by Heckmann (1975) based on the Henry Draper 
(HD) catalogue (Pickering and Cannon 1918-1924). However, Nesterov et al. (1995) 
suggested it was B9 using the extension charts of HD catalogue. The first photometric 
study, based on BVR C I C light curves, was conducted by Kleidis et al. (2008) in which, in 
addition to preliminary photometric analysis, it was suggested that HD 350731 indicates 
apsidal motion. 

In this study, the first detailed photometric and spectroscopic analysis was carried out, 
based on newly-obtained data. After descriptions of the spectroscopic and photometric 
observations, Section 3 presents spectral analysis, which contains radial velocity 
measurements of the components, preliminary orbit solution, spectroscopic light ratio, 
spectral disentangling and model atmosphere application. Analysis of multi-color light 



and radial velocity curves is given in Section 4, followed by apsidal motion analysis. The 
absolute parameters and kinematic properties of the system are detailed in Sect. 6. Finally, 
we conclude the study with the results and discussions. 


2. Observations 

New multi-color photometric observations of HD 350731 were performed at Qanakkale 
Onsekiz Mart University Ulupmar Observatory, Turkey over 11 nights in August and 
September 2012. Observations were carried out with a 60 cm Cassegrain telescope equipped 
with SB1G STL 100 IE CCD camera. One secondary minimum was also observed with a 122 
cm Cassegrain-Nasmyth telescope incorporating Apogee Alta U42 CCD camera. The data 
were collected using Bessell B , V and R c Liters. HD 350730 (A0) and HD 350727 (F5) 
were used as comparison and check star, respectively. CCD frames were reduced in the 
standard way: bias and dark frames were subtracted from the frames and then corrected for 
flat-fielding. Such reduced images were used to extract the differential magnitudes of HD 
350731. C-MunipaclJ^ code was used for these processes. Standard deviation of variations 
between the observed comparison and check stars was determined to be about 0 m .01 for all 
passbands. 

Spectroscopic observations of HD 350731 were made using the Cassegrain spectrograph 
installed on the 1.85 m Plaskett telescope at Dominion Astrophysical Observatory (DAO), 
British Columbia, Canada. The spectrograph has a spectral resolving power of about 
R = 9000. A back-illuminated SITe CCD of 1752x532 pixels (size 15 /jm) was used to 
record spectra spanning from about 4370 to 4630 A. During the observations, twenty one 
spectra for HD 350731 and several spectra for 21 Peg were taken. The spectra of the 


1 http: / / c-munipack.sourceforge.net / 
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reference star 21 Peg were used to measure the radial velocities of the components. All 
spectral data were acquired between 6 August and 6 September 2012. Spectroscopic data 
reduction was handled using the appropriate tasks of IRAlj 2 package with the following 
steps: background subtraction, division by flat field spectrum, wavelength calibration using 
Fe-Ar lamp, and normalization to the continuum. 


3. Spectroscopic Analysis 

3.1. Radial Velocities and Preliminary Orbital Solution 

In order to determine the orbital parameters of HD 350731, the radial velocities (RVs) 
of the components must be measured. Taking into account the spectral type of the system, 
given as B9 by Nesterov et al. (1995), 21 Peg (B9.5V, V r = —0.2 km s _1 ) was chosen as the 
RV standard star for both components. The cross-correlation technique (CCT) was used 
for determination of RVs of the components. Simkin (1974) reported an useful description 
of CCT which are widespread, especially for RV measurements of binary stars (e.g. Hill 
1993, Gunn et al. 1996, Frasca et al. 2002, Soydugan et al. 2007). The CCT was applied 
by FXCOR routine in the IRAF package which was developed based on the standard Tonry 
& Davis (1979) algorithm. The weights of RVs and their standard errors were calculated 
according to the usual formulas given by Topping (1972) and Tonry & Davis (1979). The 
measured RVs of both components are given in Table 1 together with their standard errors, 
which are between 4-9 km s^ 1 for the components. 

As can be seen in the light curves of the system, the secondary minimum is not 
located at the 0.5 orbital phase and indicates a displacement comparing primary minimum. 
Therefore, most probably HD 350731 has an eccentric orbit and this must be taken into 


2 http://iraf.noao.edu/ 
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Table 1: Radial velocity measurements of components of HD 350731. 


HJD 

24 50000 + 

Orbital 

Phase 

Hi 

(km s -1 ) 

h 2 

(km s _1 ) 

6173.8976 

0.0773 

- 72.9 ± 5.8 

49.1 ± 8.1 

6152.7108 

0.1200 

- 112.8 ± 5.9 

93.5 ± 7.1 

6152.7482 

0.1429 

- 126.7 ± 6.5 

111.1 ± 7.8 

6147.8861 

0.1694 

- 136.2 ± 6.2 

128.1 ± 7.0 

6152.8407 

0.1995 

- 146.9 ± 6.8 

134.8 ± 8.1 

6175.8132 

0.2488 

- 155.6 ± 6.7 

137.6 ± 8.4 

6152.9240 

0.2505 

- 152.9 ± 6.6 

137.7 ± 8.2 

6175.8962 

0.2996 

- 161.0 ± 6.8 

132.6 ± 8.1 

6180.8738 

0.3437 

- 139.1 ± 6.2 

129.4 ± 6.2 

6149.8223 

0.3536 

- 136.5 ± 6.4 

116.6 ± 8.9 

6149.9104 

0.4074 

- 97.0 ± 6.6 

89.7 ± 7.1 

6172.8880 

0.4598 

- 68.1 ± 6.5 

54.1 ± 7.4 

6177.8676 

0.5052 

- 10.9 ± 4.0 

- 

6146.9266 

0.5826 

45.7 ± 5.2 

- 62.3 ± 7.7 

6174.7848 

0.6199 

64.2 ± 6.8 

- 94.4 ± 7.2 

6151.9360 

0.6462 

99.9 ± 6.5 

- 124.1 ± 7.1 

6148.7783 

0.7151 

142.2 ± 7.7 

- 172.6 ± 7.8 

6148.8464 

0.7567 

162.4 ± 7.1 

- 181.1 ± 7.4 

6148.9094 

0.7952 

155.1 ± 6.8 

- 182.8 ± 8.5 

6171.8986 

0.8548 

136.5 ± 8.4 

- 158.5 ± 8.2 

6176.8712 

0.8958 

98.8 ± 8.8 

- 127.8 ± 6.8 
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account for the orbital solution. Using the RVs listed in Table 1 and adopted ephemeris 
from Klcidis et al. (2008), we obtained the orbital parameters of the system listed in Table 
2 in order to have input data for KOREL application and also the simultaneous analysis 
of light and radial velocity curves. The preliminary orbit parameters together with their 
uncertainties are given in Table 2. 


3.2. Spectroscopic Light Ratio 

Spectroscopic line ratios ensure us constraints on the luminosity ratios of the 
components determined from the light curve solutions at similar wavelengths, as reported 
in the study of Petrie (1939) and applied in many studies (e.g. Andersen et al. 1983, 
Southworth & Clausen 2007, Garcia et al. 2014). This also helps us to reduce the 
degeneracy, which occurs in the light curve solutions of eclipsing binaries with similar 
components for the radii of the components, if the binary are partially eclipsing system. In 
this study, we have derived spectroscopic light ratio for the system using the spectral line 
of Mgll at 4481 A. Equivalent widths (EW) of Mgll lines for the both components were 
measured using SPLOT task in the IRAF package on seven spectra taken at the phases 
out of the eclipses. The weighted mean ratio was found to be EW 2 /EW\ = 0.953T0.042, 
which can be accepted to be the light ratio of the components around 4481 A. Since the 
components of HD 350731 are almost equal and have very similar spectral properties, the 
spectral correction and wavelength dependence for light ratio are negligible. 


3.3. Spectral Disentangling 

In order to derive the atmospheric parameters of the components of HD 350731, 
individual spectrum for each component is required. The method of spectral disentangling 



Table 2: Orbital elements of HD 350731. 


Parameter Value 


To (HJD) 

Porb (day) 

V 7 (km s _1 ) 

K\ (km s _1 ) 

K 2 (km s _1 ) 
e 

co (degree) 
aisini (10 6 km) 
02 sin i (10 6 km) 
Mi sin 3 i (M 0 ) 
M 2 sin 3 i (Mg) 
q (=M 2 /M 1 ) 


2454631.4603 1 
1.635135 1 
-10.4 ± 0.7 
157.2 ± 1.3 
162.7 ± 1.3 
0.077 ± 0.007 
23.5 ± 2.3 
3.52 ± 0.04 
3.64 ± 0.04 
2.79 ± 0.05 
2.70 ± 0.05 
0.966 ± 0.015 


1 Kleidis et al. (2008) 
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of the composite spectra was invented and described by Simon & Sturm (1994). A well 
discussion and assessment on this approach was given by Pavlovski & Hensberge (2010) 
and Pavlovski & Southworth (2012). The disentangled method was also preferred in this 
study to decompose the observed spectra of the system to its components, as applied in 
some studies (e.g. Groenewegen et al. 2007, Torres et al. 2011, Lehmann et ah 2013, 
Harmanec et ah 2014). The KOREL code presented by Hadrava (1995) was used for 
spectral disentangling. Fortran code KOREL was developed based on Fourier analysis 
technique to disentangle component spectra in binary or multiple stellar systems. 

Nineteen spectra of the system obtained at DAO were used for the disentangling 
process. For the application, wavelength region 4375-4575A, which includes Mg 11 and He 
1 spectral lines, was used. The orbit properties except systemic velocity of the the system 
(Vy) could be determined with KOREL. However, the radial velocities of the components 
derived from KOREL are not independent as noted by Lehmann et ah (2013). Therefore, 
we used the orbital properties listed in Table 2 as input data for KOREL application. 
Finally, we derived the radial velocity amplitudes (K\ and K 2 ) from two different methods, 
which are compatible within better than 2 km s _1 . The fractional light contribution of 
the components determined from spectroscopic analysis was used for the application. At 
the end of the process, separate spectra of the primary and secondary components were 
obtained. Some observed composite spectra of the system at different orbital phases, 
together with the model spectra and decomposed spectrum of the components calculated 
by KOREL in the wavelength range 4450-4500 A, are shown in Figure 1. 


3.4. Model Atmosphere Application 

The physical parameters of the components of HD 350731 can be derived from analysis 


of the radial velocities of the components and light curves. On the other hand, if one 
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Fig. 1.— Fits calculated by KOREL (thick lines) to observed spectra (thin lines) at different 
orbital phases (a), and disentangled spectra of primary and secondary components (b) in 
wavelength range 4450-4500A. 
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can obtain decomposed spectra of the components, their atmospheric parameters can be 
found by line-profile fitting. The advantages of spectral disentangling method for the 
spectroscopic analysis of the eclipsing binaries’ components are presented in the studies of 
Hensberge et al. (2000) and Pavlovski & Hensberge (2005). Therefore, we obtained the 
individual spectrum of the components of HD 350731 considering these previous studies to 
make spectral analysis. The SME (Spectroscopy Made Easy) code, which enabled us to 
determine the basic atmospheric parameters of the components by matching the computed 
spectrum to the observed one was used for the analysis. SME was developed by Valenti & 
Piskunov (1996) and includes several model atmospheres to compute synthetic spectra. The 
code uses Levenberg-Marquardt algorithm to fit an observing spectrum with a synthetic 
one. The code has been applied in many studies to determine the atmospheric properties of 
stars (e.g. Torres et al. 2012, Soydugan et al. 2013). For the application, Kurucz (1993) 
model atmospheres were used and the atomic data for the spectral lines were taken from 
the Vienna Atomic Line Database (Piskunov et al. 1995; Kupka et al. 1999). 

The disentangled spectrum of both components in the wavelength region 4375-4575A 
included Mg 11 at 4481 A and He I at 4471 A spectral lines. These data were used to 
compute of the model atmospheres. The microturbulent velocities for the components 
were adopted to be 3 km s' 1 , which is an appropriate value for late-type B stars (e.g. 
Adelman 1996, Usenko et al. 2000). The surface gravities of both components were fixed 
at the dynamical values determined from solution light and radial velocity curves (see 
Table 8). After preparing SME code to fit the spectrum of the components of HD 350731, 
atmospheric parameters, namely T e ff and v sin i adjusted during the analysis have been 
determined assuming solar abundance taken from Asplund et al. (2009). The resulted 
parameters and their uncertainties estimated using the A Xmin = 1 method as described by 
Lampton et al. (1976) are listed in Table 3 . In Fig. 2, the decomposed spectra of the 
components and the synthetic spectra calculated by the best model parameters in Table 



3 are compared. The figure also shows the observed composite spectrum of HD 350731 
together with the computed spectrum, which was calculated by the model atmosphere 
parameters in Table 3, taking into account light contributions at the orbital phase of 0.2. 
The synthetic composite binary spectrum in Fig. 2 was calculated by using BinMag IDL 
visualization code (developed by O. Kochukov). As seen in Fig. 2, the synthetic spectra 
agree well with the composite binary spectrum and individual spectra of the components. 


Table 3: Model atmosphere parameters of the components of HD 350731. 



Primary 

Secondary 

Parameter 

Value 

Value 

Teff (K) 

12000 ± 250 

11830 ± 300 

logs' (cgs) 

4.25 1 

4.25 1 

v sin i (km s -1 ) 

69.2 ± 1.5 

70.1 ± 1.7 


1 Dynamical values adopted from analysis of light and radial velocity curves. 
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Fig. 2.— Disentangled spectra of components with best model atmosphere fits. Compari¬ 
son between observed composite binary spectrum and computed spectrum calculated with 
parameters in Table 3. 
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4. Analysis of Light and Radial Velocity Curves 


The BVR C light curves and radial velocities of the components of HD 350731 have been 
modeled using the version v34 of JKTEBOP code (Southworth et al. 2004, Southworth 
et al. 2005a, Southworth 2013). The code is based on the EBOP program, which was 
developed by P. Etzel (Etzel 1981, Popper & Etzel 1981). It was written in F0RTRAN77 
by J. Southwortqj and uses Levenberg-Marquardt algorithm to reach the best model. 

We preferred JKTEBOP code since it was stable, and very fast and includes various 
error-estimate algorithms. This code is also very useful for well-detached eclipsing binaries 
as tested in several studies (e.g. Ratajczak et al. 2010, Debosscher et al. 2013, Lehmann et 
al. 2013). 


HD 350731 is a partially eclipsing binary with almost identical components. In that 
case, the degeneracy may occur in determination of the radii of the components. Therefore, 
we used spectroscopic line ratio as given in the Section 3.2 to constrain the range of the 
ratio of the radii ( k ). As the first step, the luminosity ratio was adopted at 0.953 as 
determined from the spectroscopy during the analysis of B light curve since Mgll line at 
4481 A was used to derive spectroscopic light ratio. Thus, the k value was found to be 
0.993 ± 0.030. In order to check the k value and also examine the distributions of k values 
against to surface brightness ratio of the components (J 2 /J 1 ), these two parameters were 
scanned corresponding y 2 values of each solution. Then, y 2 values of the solutions were 
mapped into contours which indicate the uncertainties in Fig. 3. As shown in the figure, the 
k value of 0.993 can be seen around the lowest y 2 value. After that, the resulted k = 0.993 
was adopted during the analysis of V and R c light curves. The third light contribution 
(4) was assumed to be 0.0 after several iterations since it did not vary significantly. We 
have used the linear limb darkening coefficients as adopted parameters for the components 


3 http://www.astro.keele.ac.uk/ jkt/codes/ 
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corresponding to photometric filters used since the quadratic and logarithmic laws did 
not give better fits. The adjustable parameters are central surface brightness ratio of the 
components (J 2 /Ji), sum of the fractional radii (r 1 + r 2 ), ratio of fractional radii ( k ) for only 
B filter, orbital inclination (i), eccentricity (e sin ca, e costa), phase shift, radial velocity 
amplitudes of the components {K\ and /l 2 ) and the systemic velocity (V 7 ). 
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Fig. 3.— Contour map indicating the distribution of the results from B light curve analysis 
for the correlated parameters of surface brightness ratio (J 2 /Ji) and ratio of fractional radii 
(k) of the components. The position of the plus symbol in the figure represents minimum 
X 2 value. 
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In order to determine uncertainties for the solutions, we used task 9 of JKTEBOP 
code and calculated ler errors with Monte Carlo algorithm. The resulted parameters are 
presented in Table 4 together with their uncertainties for B,V and R c filters. The adopted 
photometric parameters, which are the weighted mean of the solutions in each filter, are 
listed in Table 5. In this table, the final orbital parameters and also their uncertainties 
calculated by JKTEBOP are also given. A comparison between observed and computed 
light curves is presented in Fig. 4 together with the residuals of the observational data from 
the best fits, while the RVs of the components together with the best fits are shown in Fig. 
5. 


5. Apsidal Motion Analysis 

The apsidal motion of HD 350731 was indicated for the first time by Klcidis et al. 
(2008). In order to establish a preliminary estimation of the apsidal motion elements based 
on the commonly-used method of O — C data analysis, we collected published minima times 
and determined two primary and two secondary minima times from our observations, as 
listed in Table 6. In total, 15 minima times were achieved. 

The analysis was made using code written by Zasche et al. (2009), which was developed 
based on the mathematical methods reported by Gimenez & Garcia-Pelayo (1983). For the 
analysis, the orbital inclination and eccentricity were adopted as i = 81°.70 from the light 
curve analysis and e=0.079 from the orbital analysis. The resulting preliminary apsidal 
motion parameters are given in Table 7, while the O — C diagram with theoretical fits is 
shown in Fig. 2. The changing rate in longitude of the periastron is Co = 0.0170 ± 0.0082 
degree cycle -1 , which corresponds to an apsidal motion period of U=92 yr. As seen in 
Table 7, the errors of the parameters are high and the parameters are not sensitive since 15 
minima times, which cover only 9 yr, were used for the analysis. However, the results can 
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Table 4: Parameters of HD 350731 obtained from analysis of multi-color light curves and 
radial velocities of components. 


Parameters 

B 

V 

R c 

n + r 2 

0.399i0.003 

0.400T0.003 

0.404i0.004 

k 

0.993T0.030 

0.993 1 

0.993 1 

n 

0.2000i0.0031 

0.2009T0.0014 

0.2028i0.0018 

?’2 

0.1986i0.0031 

0.1995T0.0014 

0.2014i0.0018 

*n 

81.83T0.09 

81.75i0.10 

81.40i0.ll 

J 2 /J 1 

0.968T0.010 

0.977i0.011 

0.976i0.015 

e sinta 

0.0349i0.0060 

0.0310i0.0020 

0.0330i0.0078 

e costa 

0.0719i0.0004 

0.0720i0.0008 

0.0717i0.0006 


1 Fixed during the analysis. 


Table 5: Final values of photometric and spectroscopic parameters for HD 350731. 


Parameters Value 


r\ 

0.2014i0.0010 

T2 

0.2000i0.0010 

in 

81.703i0.061 

L 2 /L 1 ( B) 

0.952i0.012 

l 2 /l 1 (V) 

0.961i0.011 

l 2 /l 1 (R c ) 

0.961i0.014 

e 

0.0792i0.0017 

un 

24.63i2.42 

K\ (km s _1 ) 

157.44i2.77 

K 2 (km s _1 ) 

163.14i2.22 

Vy (km s -1 ) 

-10.29il.28 
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Fig. 4.— Observed and theoretical light curves of HD 350731 in BVR C filters (a) and the 
residuals from the best fits (b). 
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Fig. 5.— RV curves of components of HD 350731 plotted versus orbital phase. Filled and 
open circles represent RVs of primary and secondary components, respectively. Solid line 
represents orbital solution for primary component and dashed line for secondary component. 
The residuals are indicated below. 






be accepted for a preliminary estimation of the apsidal motion properties and improved by 
adding new times of minima. 


6. Absolute Properties and Kinematic Behavior 

Analysis of the light curve of HD 350731 and the radial velocities of its components 
led to acquisition of its accurate absolute dimensions and physical properties for the first 
time. In order to calculate absolute parameters and their uncertainties except effective 
temperatures, which were determined from spectral analysis, we used the JKTABSDIM 
code, which is developed by Southworth et al. (2005b). In this code, the uncertainties 
are calculated with very robustly and a complete error budget is found for every output 
parameter. Derived fundamental parameters are listed in Table 8. The physical constants 
used for the calculations are given by Southworth (2011). 

The accuracy of the masses is better than 5%, while radii values of the components 
are good to about 2.5%. We calculated E(B — V ) color excess for HD 350731 from the 
dust maps of Schlafly & Finkbeiner (2011). Because the system is relatively close to 
the sun, the color excess value from Schlafly & Finkbeiner (2011) needed to be reduced 
according to the distance. The J band absolute magnitude of the system was calculated 
using the color-luminosity relation of Bilir et al. (2008). Then, the calculated color 
excess E(B — V) = 0.157 mag was reduced using the equation of Bahcall & Soneria 
(1980). The color excess value is consistent with the position of HD 350731 in the Milky 
Way. The interstellar absorption value (A v ) has been calculated to be 0.49 mag from the 
commonly-accepted formula A v = 3.1 x E(B — V). The distance of the system was found to 
be 703±34 pc based on apparent system magnitude, light ratio of the system’s components, 
and interstellar extinction values. The temperatures of the components determined from 
spectroscopic analysis are consistent with the spectral types of B8V+B8V, according to 



- 22 - 


Table 6: List of minima times of HD 350731. 


HJD 

24 50000+ 

Error 

Filter 

Epoch 

Ref. 

3594.4279 

- 

/ 

-646.5 

1 

3612.4138 

0.0001 

R c 

-635.5 

2 

3653.2919 

0.0001 

R c 

-610.5 

2 

4002.3133 

0.0001 

R 

-397.0 

1 

4299.9065 

0.0002 

V 

-215.0 

3 

4340.7838 

0.0002 

V 

-190.0 

3 

4345.6896 

0.0001 

V 

-187.0 

3 

4628.5684 

0.0002 

BVR C I C 

-14.0 

3 

4642.5482 

0.0002 

BVR C I C 

-5.5 

3 

4651.4605 

0.0004 

BVR C I C 

0.0 

3 

4665.4390 

0.0002 

BVR C I C 

8.5 

3 

6175.4058 

0.0001 

BVR C 

932.0 

this study 

6176.2986 

0.0002 

BVR C 

932.5 

this study 

6180.3108 

0.0003 

BVR C 

935.0 

this study 

6892.4851 

0.0001 

VR C 

1370.5 

this study 


x Brat et al. (2007), 2 Zejda et al. (2006), 3 Kleides et al. (2008). 
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Table 7: Apsidal motion elements of HD 350731. 


Parameter 

Value 

T 0 (HJD) 

2454651.5008 ± 0.0102 

P s (day) 

1.63511 ± 0.00001 

e 

0.079“ 

u (degree) 

3.4 ± 4.3 

cLu/dt (degree cycle -1 ) 

0.0170 ± 0.0082 

U (yr) 

92 ± 35 

P a (day) 

1.63519 ± 0.00002 


a : Adopted from orbital analysis. 
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Fig. 6.— O — C diagram for HD 350731 obtained with parameters given in Table 7. Primary 
and secondary minima times are indicated by filled and open circles, respectively. 
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calibrations given by Sung et al. (2013). 

The components of the system have similar properties since the temperature difference 
between the components was found to be AT=170 K which is smaller than la uncertainties 
for the effective temperatures of both components (see Table 3). Furthermore, the masses 
and radii values were very close, as seen in Table 8. From the photometric analysis, the 
system was found to be detached and the Roche lobe filling ratios were calculated to be 
63% and 64% for the primary and secondary components, respectively. 

In order to analyze the kinematical properties of HD 350731, we used the system’s 
center of mass velocity, distance and proper motion values. The proper motion data (/i a 
cosJ, /i< 5 ) = (1.0±0.7, -8.7±0.8) mas yr -1 were taken from the Fourth US Naval Observatory 
CCD Astrograph Catalog (UCAC4; Zacharias et al. 2013), whereas the center of mass 
velocity 1% = —10.1 ± 0.4 km s -1 and distance d=703±34 pc were obtained in this study. 
The system’s space velocity was calculated using Johnson & Soderblom’s (1987) algorithm. 
To obtain the space velocity precisely, first-order galactic differential rotation correction 
was taken into account (Mihalas & Binney 1981). The differential rotation corrections were 
calculated as 16.23 and 0.89 km s _1 and applied to U and V space velocity components, 
respectively. The W velocity is not affected in this first-order approximation. As for 
the local standard of rest correction (LSR), Co§kunoglu et al.’s (2011) values (U, V, 

ID) = (8.50, 13.38, 6.49) km s -1 were used and the total space velocity of HD 350731 was 
obtained as S tot = 15.30±3.84 km s -1 . The corrected space velocity components are ( U , 
V, 1 D ) = (6.10±2.39, -9.26±1.56, -10.54i2.57) km s -1 . The total space velocity and space 
velocity component values are in agreement with young-disc stars (Leggett, 1992). 

To determine the population type of HD 350731, we used Dinescu, Girardi & van 
Altena’s (1999) N-body code and obtained the galactic orbit of the system. In this code, 
the timescale in generating the orbits was assumed to be 3 Gyr, and the calculation steps 



Table 8 : Astrophysical properties of HD 350731. 


Parameter 

Primary 

Secondary 

Mass (M 0 ) 

2.91±0.13 

2.80i0.14 

Radius (R e ) 

2.11T0.05 

2.07i0.05 

Temperature (K) 

12000T250 

11830i300 

log L (L 0 ) 

1.92T0.04 

1.88i0.05 

logs- (cgs) 

4.25T0.02 

4.25i0.02 

Orbital period (day) 

1.63511±lxl0" 5 

Orbital separation (R 0 ) 

10.43i0.02 

Mass ratio 

0.965±0.004 

Systemic velocity (kms' 1 ) 

-10.1i0.4 

Distance (pc) 

703T34 

V (mag) 

9.60 a 

M B oi (mag) 

-0.05i0.ll 

0.05i0.12 

BC (mag) 

-0.66 fc 

-0.62 6 

M v (mag) 

0.61i0.11 

0.67i0.12 

Measured v sin i (kms^ 1 ) 

69.2il.5 

70.1il.7 

Pseudo-synchronous v sin i (kms“ 

- 1 ) 67.5i0.3 

67.0i0.3 

Synchronous v sin i (kins' 1 ) 

65.4i0.3 

64.1i0.3 

Age (Myr) 

120i35 

“:SIMBAD Database, b : 

Sung et al. (2013) 
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were 2 Myr. The 3 Gyr timescale was assumed so that precise orbits would be created, 
even though this is longer than the nuclear time scale of early-type stars. The orbits of HD 
350731 on the X — Y and X — Z planes around the galactic center are shown in Fig. 7. 
The system’s apogalactic (R max ) and perigalactic (. R m in ) distances obtained were 7.72 and 
7.09 kpc, respectively. According to N-body code, the maximum vertical separation from 
the galactic plane of the system is | z max | = 100 pc. The following formulae were used to 
derive the planar (e p ) and vertical (e v ) eccentricities: 


e 


p 


R 

R 


max 


max 


Rr 


+ Rrt 


( 1 ) 


c v 


(| %max | | ^min |) 

Rrn 


( 2 ) 


where R m is the mean of R m i n and R max ■ The planar and vertical eccentricities were 
calculated as e p = 0.04 and e v = 0.01, respectively. These eccentricities show that HD 
350731 is in a circular orbit around the mass center of the Galaxy and that it belongs to 
the young thin-disc population. 


7. Discussion and Conclusion 

Double-line detached eclipsing binaries (DBs) are valuable sources for determining the 
precise fundamental properties (mainly masses and radii) of stars. A recent catalogue of 
this type of binaries was published by Eker et al. (2014). It consists of 257 DBs; their 
388 component stars have better than 5% accuracy in their masses and radii. When one 
examines the basic properties of the DBs in the catalogue (especially mass ratio, mass, 
radius and temperatures of components), it can be seen that there is no detached eclipsing 
binary system which has absolute properties similar to HD 350731. The spectral type 
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Fig. 7.— Orbital motion of HD 350731 on projections of X — Y (a) and X — Z (b) planes 
around galactic center for 3 Gyrs. 

















of the system was determined to be B8V+B8V, while the masses were derived to be 
Mi = 2.91 ± 0.13 M 0 and M 2 = 2.80 ± 0.14 M 0 for the primary and secondary components, 
respectively. Therefore, DBs with mass ratio in the range of 0.95-1.0 and spectral type O 
and/or B need to be further studied in order to fill the gap in this parameter range. This is 
necessary to understand the evolution and structure of early-type stars. 

The mass-luminosity relation (MLR) was updated recently by Eker et al. (2015) using 
DBs data. In this study, they identified four mass domains and derived MLRs for these 
mass ranges. For comparison, using the MLR given in the mass range of 2.4 M 0 and 7M 0 , 
we calculated luminosity values on a logarithmic scale to be 1.96 L 0 and 1.89 L 0 , for the 
primary and secondary components, respectively. The predicted values from the updated 
classical MLR agree with the values derived from Stefan-Boltzmann law (see Table 8). 

Spectral data enables us to find temperatures and projected rotational velocities of 
the components. For the model atmosphere application, the disentangled spectra of the 
components obtained from KOREL analysis were used. As a results the temperature of the 
primary and secondary components was found to be 12000 K and 11830 K, respectively. 
The surface gravity values (logg li2 ) calculated from the masses and radii of the components 
were adopted during spectral analysis. The projected rotational velocities (vsini) of 
the components were measured to be 69.2 kms -1 and 70.1 kms -1 for the primary and 
secondary components, respectively. Within errors the measured v sin i values agree with 
the pseudo-synchronous velocities of the components in Table 8, which were calculated on 
the basis of the formulations by Hut (1981). 

The orbit of the system is slightly eccentric determined from the orbital solution and 
also an analysis of the multi-color light curves and radial velocities of the components. The 
secondary minimum can be seen to have shifted from the orbital phase of 0.5. This is an 
indication of apsidal motion in the system. We collected 15 minima times together with 
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newly-measured eclipse times to study the apsidal motion by means of O — C analysis. The 
apsidal motion parameters could not be determined accurately since the observed eclipse 
times were not-well covered. The apsidal motion was found at a rate of uj — 0.0170 ± 0.0082 
degree cycle -1 , which corresponds to an apsidal motion period of about 92 yr. 

The locations of the components of HD 350731 are seen in the Hertzsprung-Russell 
(HR) diagram and plane of M — R in Fig. 8. Zero Age Main Sequence (ZAMS) and 
the evolutionary tracks for the exact masses of the components, and isochrones for solar 
chemical composition, are taken from the Yonsei-Yale (Y 2 ) series of Yi et al. (2001). As 
seen in Fig. 8, the components are a bit away from the ZAMS line and the age of the 
system was estimated to be 120±35 Myr from Y 2 isochrones. Kinematic analysis indicated 
that the system leaves the galactic plane just about 100 pc during its movement on galactic 
orbit. This is evidence for membership of the thin-disc population by HD 350731. The 
positions of the components in the HR diagrams are in agreement with the evolutionary 
tracks for the masses of 2.91 M e and 2.80 M 0 . 

As a general conclusion, we can mention that the photometric and spectroscopic 
analysis of HD 350731 leads us to extend the database of DBs with intermediate mass 
components having similar properties. This is important since in the catalogue of DBs 
by Eker et al. (2014) there are only a few systems with a mass ratio close to q — 1.0 
and components having masses greater than 2.8 M 0 (e.g. rj Mus, V799 Cas, V906 Sco). 
The unevolved binary system HD 350731 belongs to the young thin-disc population of the 
Galaxy with an estimated age of 120±35 Myr. In order to make abundance analysis in 
detail, and also verify population type of the system, more high resolution spectroscopic 
data is needed. New eclipse times are also required to enlarge the time interval for better 
analysis and to confirm the apsidal motion parameters. 

This research was supported by the Scientific and Technological Research Council of 
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log T eff (K) 



log (M/M @ ) 


Fig. 8.— Locations of primary and secondary components of HD 350731 in log L - log T e ff 
(a) and log M -log R planes (b). Evolutionary tracks for masses of 2.91 M 0 and 2.80 M 0 (solid 
lines), isochrones for ages of 80 Myr, 120 Myr and 150 Myr (dashed lines), and ZAMS for 
solar chemical composition adopted from Yi et al. (2001). 
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